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1.  Introduction  and  Abstract 


By  a  cell  we  mean  either  a  nonempty  closed  polyhedral  convex 
set  or  a  nonempty  closed  solid  ball.  Our  concern  is  with  solving  as 

t 

linear  or  convex  quadratic  programs  special  cases  of  the  following 
two  problems. 

Optimal  Containment  Problem  (OCP): 

Let  St  be  a  finite  union  of  cells  and  a  finite  intersection 
of  cells.  Find  the  smallest  positive  scale  sty  of  for  which 

some  translate  s<&  +  t  contains  St. 

inf  iraum :  s 
8,t 

subject  to:  +  t  ^  St,  s  >  0  . 

Optimal  Meet  Problem  (OMP): 

Let  St  and  3^  for  q  =  1,  ...  ,  p,  be  each  an  intersection 
of  cells.  Find  the  smallest  positive  scale  sSt  of  St  for  which 
some  translate  sSt  +  t  meets  every  &  for  q  =  1,  ...  ,  p. 

s 

{sSt  t)  O  3^  ^  ,  q  —  1| .  * .  jP  s  >  0  . 

The  general  OCP  and  OMP  are  well  beyond  our  reach  but  serve  as 
useful  overviews.  Depending  upon  the  compositon  of  the  St's  and 
^»3  as  unions  and  intersections  of  cells  and  the  representation  of 
the  cells  we  can  or  cannot  formulate  the  problems  as  linear  or  convex 
quadratic  programs. 


OMP  < 


infimum: 


a,t 


subject  to: 


Our  initial  interest  in  OCP  and  OMP  originated  with  "engineering 
design"  through  interior  solution  concepts  for  convex  sets,  see  van  der 
Vet  [ 9  ]  and  Director  and  Hachtel  [  1  ];  also  see  Eaves  and  Freund  [  2  ]. 

2.  Preliminaries 

Most  of  the  notation  we  use  is  standard.  Let  IRn  be  n-dimen- 
sional  Euclidean  space.  By  ||*||,  we  mean  the  Euclidean  norm,  x  •  y 
and  x  o  y  represent  inner  and  outer  product,  respectively.  Let 
e  =  (l,  1,  ...,  l)  where  its  length  is  dictated  by  context.  We  define 
inf  <(>  =  +“  as  usual,  but  sup  $  ^  0  for  the  purposes  of  our  presen¬ 
tation.  By  a  convex  program  we  mean  a  program  of  the  form 

f(x) 

6^(x)  ^  0  >  i  s  1,  •••»  B  , 

where  all  f,  g^  are  convex  functions,  and  m  is  finite.  If,  in 
addition,  each  g^  is  affine,  f(x)  *  xQx  +  q  •  x,  and  Q  is 
positive  semi-definite,  then  we  call  P  a  quadratic  program.  Further¬ 
more,  if  Q  is  zero,  then  we  call  P  a  linear  program. 

Let  g  be  a  set  in  lRn.  We  denote  by  tng(^)  the  smallest 
vector  subspace  of  3Rn  for  which  some  translate  contains  g .  We 
denote  by  rec (?)  the  recession  set  of  £  ,  that  is,  the  set 

(z  c  B0! 3x  e?  such  that  x  +  crz  e  £  for  any  a  >  0}  . 


p 


r  minimum: 
x 


subject  to: 


We  also  make  use  of  the  following  variation  of  Far  leas'  Lemma 


Lemma:  Suppose  that  the  system  of  inequalities 

Ax  <  b 

has  a  solution  and  that  every  solution  satisfies  cx  <  d.  Then  there 
exists  X  >  0  such  that  XA  =  c  and  Xb  <  d.  0 

The  manner  in  which  cells  are  represented  is  crucial  to  our 
formulations.  We  assume  all  cells  are  in  3Rn.  We  define  an  H-cell 
to  be  a  cell  of  the  form 

(x|Ax  <  b} 

as  it  is  represented  by  half-spaces;  it  is  assumed  (A,b)  is  given. 

A  cell  of  the  form 

{xjx  *  UX  +  Vn  :  eX  *  1,  X  >  0,  n  >  0} 

is  defined  to  be  a  W-cell  as  it  is  a  weighting  of  points;  it  is  assumed 
(U,V)  is  given.  Of  course,  every  H-cell  can  be  represented  as  a  W- 
cell  and  vice  versa.  However,  we  shall  suppose,  and  typically  rightly 
so,  that  the  computational  burden  of  the  conversion  is  prohibitive, 
see  Mattheiss  and  Rubin  [8]  for  H  to  W.  Thus  we  shaUL  regard  H- 
and  W-cells  as  quite  distinct.  A  B-cell  is  defined  to  be  the  ball 

(x|  ||c  -  x||  <  r}  ; 

it  is  assumed  that  the  center  c  and  radius  r  >  0  are  given. 


Let  J  be  a  cell,  and  let  s  >  0  be  a  scale  of  g  and 
t  a  translate.  If  g  is  an  H-cell,  (x|ax  <  b},  then 
s^+  t  =  (x|Ax  <  bs  +  At}.  If  g  is  a  B-cell,  [x|  ||c  -  x[|  <  r}, 
a  g*  +  t  =  (x|  || (sc  +  t)  -  x||  <  sr}.  And  if  %  is  a  W-cell, 

(x|x  =  UX  +  Vji,  eX  ■  1,  X  >  0,  (i  >  0},  then  s^  +  t 
=  (x|x  *  (sU  +  t®  e)x  +  V| i,  eX  =  1,  X  >  0,  n  >  0}  or  equivalently 
aj  +  t=  (x|x=t+UX  +  Vji,  eX  =  s,  X  >  0,  |i  >  0). 

(7 

To  describe  special  cases  of  OCP  and  CM*  we  shall  use  notation 
as,  for  example,  (HB1,  WB_)  which  denotes  that  0C  is  composed  of 
any  finite  number  of  H-cells  and  one  B-cell,  and  that  or  each 

are  composed  of  any  number  of  W-cells  and  B-cells  but  all  B-cells 
have  the  same  radius.  If  B  is  not  subscripted  by  "="  or  "1”, 
then  any  finite  number  may  be  employed  and  the  radii  may  vary.  Thus, 
again,  for  example,  (HWB,  HWB)  describes  the  most  general  case  of 
OMP  or  OCP. 

Consider  the  following  three  programs 

s 

1^  -  (sc  +  t ) [|  +  d  <  sr,  i  =  1,  ...,  m  . 
f 

l|w±  -  x||2  <  f ,  i  -  1,  . . . ,  m  . 
x  *  x  -  a 

Wi  *  wi  ”  2wi  '  x  +  a  1  i  “  1,  ...,  ra  . 


(<U)  { 


f  v1  =  infimum: 
s,  t 

subject  to: 


(02) 


v_  =  infimum: 
"  f,x 

subject  to: 


(03) 


v,  *  infimum: 
5  a,x 

subject  to: 


6 


1 


Assuming  r  is  non-zero ,  we  show  that  solving  any  one  yields 
solutions  to  the  other  two.  Let  V7  denote  nonnegative  square  root. 

Equivalence  of  (Ql)  and  (02);  If  (s,t)  and  (s,t)  are 
feasible  for  (Ql)  with  s  <  s  (respectively:  s  <  s),  then 

(f,x)  =  ((sr  +  d)^,  sc  +  t)  and  (f,x)  =  ((sr  +  d)^,  sc  +  t)  are 

feasible  for  (Q2)  with  f<f  (respectively:  f  <  f).  If  (f,x) 

and  (f>x)  are  feasible  for  (Q2)  with  f  <  f  (respectively: 

f  <  f),  then  (7,t)  =  ((Y?  +  d)/r,  x  -  sc)  and  (s,t)  = 

((V?  +  dVr>  x  -  sc)  are  feasible  for  (Ql)  with  s  <  s  (respectively: 
s  <  s). 


Equivalence  of  (Q2)  and  (Q3) :  If  (f,x)  and  (f,x)  are  feasible 
for  (Q2)  with  f<f  (respectively:  f<f),  then  (a,x)  (x*x-f,  x) 
and  (a,x)  =  (x  •  x  -  f,  x)  are  feasible  for  (Q5)  with  x  •  x  -  a 

=f<x*x-a=f  (respectively:  x*x-a  =  f<x»x-a  =  f).  ' 

If  (a#:D  and  (a,x)  are  feasible  for  (Q5)  with  x*x  -  a  <  x*x  -  a 

(respectively:  x*x-a<x*x-a),  then  (f,x)  =  (x  •  x  -  a,  x) 

and  (f,x)  =  (x  •  x  -  a,  x)  are  feasible  for  (Q2)  with  f  »  x  •  x  -  a 

<  f  =  x  •  x  -  a  (respectively:  f  =  x*x-a<f  =  x*  x-a). 

We  thus  have  the  following  result 
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Lemma  2.1  (Equivalence  of  (Ql)  and  (Q3)):  For  r  non-zero 

(i)  If  (s,t)  is  feasible  or  optimal  for  (Ql),  then  (a,x)  = 

((sc  +  t)  •  (sc  +  t)  -  (sr  -  d)^,  sc  +  t)  is  feasible 

or  optimal  for  (03),  respectively. 

(ii)  If  (a,x)  is  feasible  or  optimal  for  (03),  then  (s,t)  = 

((d  +  1/x  •  x  -  a)/r,  x  -  (d  +  -^x  •  x  -  a)c/r)  is 

feasible  or  optimal  for  (Ql),  respectively.  0 


Consequently  (Ql)  can  be  solved  via  the  quadratic  program 
(03) •  Note  that  (03) ,  and  hence  (Ql) ,  always  has  a  unique  optimal 
solution. 


3*  The  Optimal  Containment  Problem  (OCP) 

Let  be  a  finite  union  of  cells  and  ‘3/  be  a  finite  inter¬ 
section  of  cells.  The  optimal  containment  problem  can  be  written  as: 


or  as 


CCP1 


{ 


z.  =  supremum: 
s,t 

subject  to: 


s 

sSC  +  t  C 


s  >  0 


zp  =  infimum: 
0CP2  l  *  s ,t 


subject  to:  3C  C  s&  +  t,  s  >  0 


8 


ju*. 


1 1  atffeV  f  •  ifllffiftr--'  • 


1 1-  jM  ii  ir 


Our  first  results  concern  the  equivalence  of  0CP1  and  0CP2. 


Lemma  3.1  (Equivalence  of  0CP1  and  0CP2) : 

(i)  (Solutions)  (s,t)  is  a  feasible  or  optimal  solution  to 

0CP1  if  and  only  if  (l/s,  -t/s)  is  a  feasible  or  optimal 
solution  to  0CP2. 

j 

(ii)  (Feasibility)  The  following  are  equivalent :  ; 

! 

(a)  0CP1  is  feasible  ' 

(b)  0CP2  is  feasible  - 1 

(c)  tng(^)  C  tng(S^)  and  rec(^)  C  rec(t^). 

(iii)  (Attainment)  The  following  are  equivalent: 

(a)  0  <  <  +  «>  (c)  0CP1  has  an  optimum 

(b)  0  <  Zg  <  +  »  (d)  0CP2  has  an  optimum. 

(iv)  (Non-attainment)  The  following  are  equivalent: 

(a)  z1  =  +  <» 

(b)  z2  =  0 

(c)  SC  +  t  C  rec(S/')  for  some  translate  t. 

For  a  specific  realization  of  0CP1  or  0CP2,  SC  and 
will  be  given  in  the  forms 
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St' =  t  U  H.]  U  [  U  W.]u[u  P.] 
hea  ieP  1  je7  2 


where 


\  >  0, 


whose 


a  =  $ 
H-,  W- 

Case 


0CP2. 


V  =  [  n  il  ]  n  [  n  w  ]  n  [  n  Bm]  , 

kep  ieff  mex 


H(.)  =  klA(.)x  <  b^},  =  (x|x  =  U^\  +  V^p,  e\  =  1, 

H  >  0},  and  ^  =  (x|  llc^.j  -  x||  <  r^j).  For  a  given  set 

:  tx|A(.)X  <b(>)},  we  define  a(.)  to  be  the  column  vector 

til  tli 

q  component  is  the  (Euclidean)  norm  of  the  qX  row  of  A^j, 

We  begin  with  case  (HWB,  H)  of  OCP  which  corresponds  to 

and  t  —  p,  that  is,  SC  is  a  union  of  any  finite  number  of 

,  and  B-  cells  and  is  an  intersection  of  H-cells. 


(HWB,  H)  of  OCP  is  a  linear  program 

We  treat  the  optimal  containment  problem  (HWB,  H)  through 
The  formulation  as  a  linear  program  is 


=  minimum:  s 

s ,  t ,  A 

subject  to:  A^A^  =  Ak 

hea. 

kep 

Ahkih  S  V  +  V 

hea. 

kep 

Vi'-V'V 

i  £  P, 

kep 

Vh’s0 

i  €  P, 

kep 

Vj +  Vj  i  V  +  V 

j  €  V, 

kep 

hea. 

kep 

s  >  0  . 

10 


is 


over  a  decade. 


Case  (Wt  HW)  of  OCP  is  a  linear  program 
Formulated  through  0CP1,  we  have: 


=  maximum  :  s 


subject  to:  sU^  +  toe  =  U^A^  +  V^ir^ 


Vi  =  Vi! 


i  e  p,  l  e  a 
i  e  p,  l  e  a 
i  e  p,  k  e  p 
i  e  p,  k  e  p 


Aj^sl^  +  toe)  <  bkoe 

W  £° 

eAii  -  e>  4i/  >  o-  "i  i  >  o,  nu  >  o  1  =M.. 

s  >  0  . 

Case  (B^,  B_)  of  OCP  is  a  quadratic  program 

let  (c,r)  be  the  given  center  and  radius  of  the  ball  SC. 
Treating  the  optimal  containment  problem  through  0CP1,  we  have: 

z1  =  maximum:  s 
1  sf  t 

subject  to:  ||c  -  (sc  +  t)||  +  sr  <  r  n  e  t 

ra  —  n 

s  >  0  . 
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- . ~  •  -  - . 


If  j  i.e*>  the  intersection  of  the  Bffl  is  not  empty,  then 
the  constraint  s  >  0  is  superfluous,  and  can  be  dropped.  Since  all 
rn  are  equal,  the  above  program  is  seen  to  be  an  instance  of  (Ql)  and 
hence  can  be  solved  via  the  quadratic  program  (Q3).  Note  that  if 
the  optimal  solution  to  the  program  (Q3)  returns  a  negative  value 
of  s,  then  3^  =  $.  A  variation  of  this  problem  was  first  shown 
to  be  a  quadratic  program  by  Gale  [  4  ]. 

Case  (B_,  B^)  of  OOP  is  a  quadratic  program 

Here  we  let  (c,r)  be  the  given  center  and  radius  of  the  ball 
“0/ .  Formulated  through  0CP2,  the  optimal  containment  problem  is 
written  as 

z_  =  minimum:  s 
2  s,t 

subject  to:  Hc^  -  (sc  +  t)||  +  r..  <  sr  j  e  7 

s  >  0  . 

Note  that  the  constraint  s  >  0  is  superfluous,  and  can  be  omitted. 

As  this  program  is  a  realization  of  (Ql),  it  is  solvable  as  the 
quadratic  program  (Q3)  for  r  >  0. 

The  special  case  of  (B_,B^)  where  all  r^  =  0  is  the  problem 
of  finding  the  smallest  ball  covering  the  points  c^,  j  €7  and  has 
been  treated  by  Elzinga  and  Hearn  [J]  and  Kuhn  [7]. 


The  case  (W,  B^)  of  OCP  is  a  quadratic  program 

(W,  B1)  of  OCP  is  a  special  case  of  (B_,  B^)  just  discussed 
when  rec(&')  =  )  (else  OCP  is  infeasible),  since  all  =  0  and. 

each  column  of  the  can  be  considered  as  the  center  of  a  B-cell 

with  radius  zero.  Thus  (W,  B^  of  OCP  is  solvable  through  the 
quadratic  program  (Q3). 

Other  cases  of  OCP 

Cases  (WB,  HB)  and  (W,  HWB)  of  OCP  can  be  formulated  as 
convex  programs  using  the  logic  already  employed;  however,  we  have  been 
unable  to  formulate  either  case  as  a  quadratic  or  linear  program.  As 
regards  all  other  cases  of  OCP,  we  are  convinced  that  their  formu¬ 
lation  as  a  convex  program,  much  less  a  quadratic  or  linear  program, 
cannot  be  accomplished.  The  reason  for  this  is  that  the  problem  of 
testing  3C  C  aj ,  where  either  (i)  &'  is  an  H-cell  and  gf  Is  a 

W-cell,  (ii)  is  an  H-cell  and  is  a  B-cell,  or  (iii)  9t'  is 

a  B-cell  and  ^  is  a  W-cell,  appears  to  be  intractable  without 
conversion  of  the  polyhedra  from  H-cell  to  W-cell  or  vice  versa. 

4.  The  Optimal  Meet  Problem  (OMP) 

Let  9'  and  q  *»  1,  ...,  p,  each  be  a  finite  intersection 

of  cells.  The  optimal  meet  problem  can  be  written  as: 

f*  v1  *  infimum:  s 
1  s,t 

subject  to:  (s^T  +  t)  H  ^  ^  <*> 


0MP1 


q  =  1,  . . . ,  p  , 


QMP2 


i 


"  v9  -  supremum:  s 
s,t 

<  subject  to:  3l  n  (s  +  t)  ^  ♦  q  =  1,  ...»  p 
^  s  >  0  . 

The  following  result  concern  the  equivalence  of  0MP1  and  0MP2. 

Lemma  4.1  (Equivalence  of  0MP1  and  0MP2) : 

(i)  (Solutions)  (s,t)  is  a  feasible  or  optimal  solution  to 

0MP1  if  and  only  if  (l/a<  -t/s)  is  a  feasible  or  optimal 
solution  to  0MP2. 

(ii)  (Feasibility)  The  following  are  equivalent: 

(a)  0MP1  is  feasible 

(b)  0MP2  is  feasible 

(c)  For  some  translate  t,  (t  +  tng (£'))  O  y  f 

q  *  i»  •••!?• 

(iii)  (Attainment)  The  following  are  equivalent: 

(a)  0  <  <  ■  (c)  0MP1  has  an  optimum 

(b)  0  <  Vg  <  «  (d)  0MP2  has  an  optimum. 

(iv)  (Non-attainment)  The  following  are  equivalent: 

(a)  v'1  ■  0 

(b)  v2  »  +  » 

(c)  For  some  translate  t,  (rec (.^T)  +  t)  n 

q  *  •••>  P  • 
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For  a  specific  realization  of  0MP1  or  0MP2,  and  ^ 

will  be  given  in  the  form: 

JT-  [  n  h  ]  o  {  n  w  ]  n  [  n  bJ 

hea  n  teg  1  Jcr  13 

d  -  [  n  h.  ]  n  [  n  if  i  n  [  n  b  ]  q  *  l,  . . . ,  p  , 

VP1  q  V°q  9  V\  9 

where  ^(-)  an<*  ®(.)  &re  as  ^  section  3.  Cfctr  solvable  cases 

are  as  follows. 


Case  (HW,  HW)  of  OMP  is  a  linear  program 

Treated  through  0MP1,  case  (KW,  HW)  is  the  linear  program: 

v,.  »  minimum  :  s 
s  f t,X  jX;4 


subject  to: 


Vq  ±  bh“  +  V 


V 


-  \ 
q  q 


Xq“  t+Ui\q+ViUiq 


*q  *  Ui  N  +  V/ 

q  q  q  q 


e\q  “  »'  No  i  °»  No  ^  0 


iq  ~  ’  iq 


e\.  *  1,  \£  >  0,  \Lt  >  0 

q  q  q 


h  €  Cl*  q  * 

kq  e  pq'  9  * 

1  e  fit  q  “  1/*»*»P 

q  *  1>«»**P 

l  e  fit  q  ■ 

eq  €  *V  9  *  1'***,p 


8  >  0  . 
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Case  (B^,  B^_)  of  OMP  is  a  quadratic  program 

This  is  the  case  where  each  <&  is  a  ball  of  given  center  c^ 

and  radius  r  with  r,  =  •••  =  r  >0.  Let  (X  be  the  ball  with 

q  1  q  - 

center  c  and  positive  radius  r,  and  we  proceed  through  0MP1. 

The  formulation  is 


v  =  minimum: 
1  s,t 


subject  to: 


(sc  +  t)||  <  sr  +  r^  , 
q  =  1,  ...,  p 


s  >  0  . 


P 

If  H  =  $  (otherwise  OMP  does  not  attain  its  minimum),  then 
q=l  q 

the  constraint  s  >  0  is  superfluous  and  can  be  emitted.  Furthermore, 
since  all  r^  are  equal,  the  above  program  is  an  instance  of  (Ql) 
and  hence  can  be  solved  via  the  quadratic  program  (<13) . 


Other  case 3  of  OMP 

Rote  that  the  most  general  case  of  CM1,  namely  (HWB,  HWB), 
can  be  formulated  as  a  convex  program  using  the  logic  employed  herein. 
However,  we  have  been  unable  to  formulate  any  case  of  OMP  other  than 
the  above  two  cases  as  a  quadratic  or  linear  program. 
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5.  Remark 

Our  final  remark  concerns  the  interrelated  issues  of  computational 
complexity,  the  conversion  of  H-  to  W-cells  and  vice  versa,  and  our 
division  of  problems  solvable  as  quadratic  or  linear  programs  from 
other  convex  programs.  In  (53  and  [6],  it  is  shown  that  linear  and 
(convex)  quadratic  programs  are  solvable  in  polynomial  time.  The 
conversion  of  an  H-  to  W-cell,  or  vice  versa,  is  an  exponential  problem. 
To  see  this,  consider  the  sets  A  (x  €  Rn|  ||x||«»  <  l)  and 
Q)  A  {x  €  Rn|  HxL  <  l).  V  ,  as  an  H-cell,  can  be  represented  by 
2n  half spaces,  but  as  a  W-cell  it  requires  the  enumeration  of  at 
least  2n  (extreme)  points.  Q  ,  as  a  W-cell,  can  be  represented 
by  2n  (extreme)  points,  but  as  an  H-cell  it  requires  the  enumeration 
of  at  least  2n  halfspaces.  We  thus  see  that  our  distinction  of 
H-cells  and  W-cells  as  different  entities  is  consistent  from  the 
standpoint  of  the  solvability  of  the  quadratic  and  linear  programs 
contained  herein. 
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